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Abstract 

We study synchronization of oscillators that are indirectly coupled 
through their interaction with an environment. We give criteria for 
the stability or instability of a synchronized oscillation. Using these 
criteria we investigate synchronization of systems of oscillators which 
are weakly coupled, in the sense that the influence of the oscillators on 
the environment is weak. We prove that arbitrarily weak coupling will 
synchronize the oscillators, provided that this coupling is of the 'right' 
sign. We illustrate our general results by applications to a model of 
coupled GnRH neuron oscillators proposed by Khadra and Li [14], and 
to indirectly weakly-coupled A — uj oscillators. 

PACS: 05.45.Xt; 87.18. Gh; 87.18.Hf 
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1 Introduction 



The aim of this work is to investigate the dynamics of systems of oscillators 
which are globally coupled through an environment. An important example 
of such systems is populations of cells in which oscillatory reactions are taking 
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place [H QUI [151 ETJ E5] , which 'communicate' via chemicals that diffuse in the 
surrounding medium. The ability of thousands of cells to synchronize their 
periodic activity is crucial for the generation of macroscopic oscillations, like 
circadian periodicities [2]. 

Consider a system of n identical dynamical systems ('oscillators'), described 
by the differential equations 



where Xk G M d is the state of the k-th oscillator, y G M p is the state of the 
environment, and / : R^xK^ R d is smooth. The dynamics of each oscil- 
lator thus depends on the state of the environment. An additional equation 
describes the dynamics of the environmental variable y 



where the smooth function g : M p — > MP represents the intrinsic dynamics of 
the environment, and the smooth function h : M d x W p — > MP represents the 
effect of the oscillators on the environment. The state of the oscillators thus 
influences the dynamics of the environment. 

In the case of biological cells, Xk would be a vector whose components are 
the concentrations (in moles per unit volume) of various biochemical species 
in cell k, and y a vector of concentrations of various biochemical species in 
the exterior of the cells. The parameter (3 is the ratio of the total intracel- 
lular volume to the volume of the environment: if V ce u is the volume of an 
individual cell and V ext is the volume of the external environment, 



In the case of biological cells, the interaction of cells and environment may 
occur through the diffusion and transport of chemical species across the cell 
membranes, and through the effects of the activation of receptors on the 
cell membrane. A variety of modeling studies of biochemical systems of 
oscillators coupled through an environment, described by equations of the 



form(0),([L2]), can be found in [U El H EJ 1 El UM TH E21 ESI ES E3 [28]. The 



framework presented above thus unifies many models of particular systems, 



x k = f(x k ,y), 1 < k < n, 



(1.1) 




(1.2) 
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and allows us to obtain some basic analytical results which apply to all of 
them. 

Since the state of each of the oscillators influences the environment, and 
the state of the environment in turn influences the oscillators, we obtain an 
indirect coupling of the oscillators. We are interested in studying the capacity 
of this indirect coupling to induce synchronization of the oscillators. When 
this happens, in the biochemical context, one refers to the relevant species 
which diffuse in the environment as 'synchronizing agents'. 

By a synchronized oscillation of the system (11.11) . (11.21) we mean a periodic 
solution with all the Xk identical 

xi(t) = x 2 (t) = ... = x n (t) = x(t). (1.4) 

The periodicity means that there exists a T > such that 

x(t + T)=x(t), y(t + T)=y(t), Vi G R. (1.5) 

Therefore, substituting (jl.4p into (II. ip . (11.21) . we see that a synchronized os- 
cillation corresponds to a periodic solution of the system 

x = f(x,y), (1.6) 

y = g{y)+/3h{x,y). (1.7) 

The original system (jl.ip . (jl.2p is nd + p dimensional, whereas the system 
(11.6p . (jl.7j) is only d + p dimensional. Since (jl.6p . (jl.7j) do not depend on n, 
a periodic solution of (ll.6p . (ll.7p gives rise to a synchronized oscillation of 
(rrrp .fTO] for anyn> 1. Note that system (JUSP , (HZZD is simply (TO) . (TO) 
for the case n — 1, so it describes the behavior of a single oscillator placed 
in the environment, and we can therefore refer to it as the single-oscillator 
system. 

In order for a synchronized oscillation to be observable, it must be stable 
in the sense that it is asymptotically approached starting from an open set 
of initial conditions (the precise definition is recalled in section [2]). A cru- 
cial point must be made here: the stability of a synchronized oscillation 
x\(t) = ... = x n {t) = x(t), y{t) refers to its stability as a solution of the full 
system (II. II) . (II. 21) . and does noHollow from the stability of the corresponding 
solution (x(t),y(t)) as a solution of the single-oscillator system (ll.6p . (ll.7p . 
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When we say that synchronization occurs, we mean that there exists a syn- 
chronized oscillation which is stable as a solution of (jl.lj) . (11.20 . A criterion 
for the (in)stability of a synchronized oscillation will be derived in section [3j 

In order to illustrate some of the collective phenomena observed in systems 
of the form fll.ip . fll.2p we display, in figured! the results of numerical simu- 
lations, in which we took three Van der Pol oscillators coupled through an 
environment. The dynamics of the oscillators are defined by 

Vk + {v 2 k - l)vk + Vk = y, 1 < k < 3, 

which we can convert to a first-order system fll.ll) . with d — 2, p — 1, n — 3, 
where x = (v,w) and 

f(x, y) = O, (1 - v 2 )w -v + y). 

We assume now that the environmental variable evolves according to fll.2p . 
with 

g{y) = -y, K x ^v) = v - 

The system of the three oscillators coupled through the environment is thus 
described by the equations 

Vk + (vl - l)v k + v k = y, 1< k<3 

y = -y+ + v 2 + v 3 ). 

In the top part of figure [H we set (3 = —0.5, and plotted the values Vk(t) 
(1 < k < 3), starting at some arbitrary initial conditions. One sees that the 
three oscillators synchronize, so that within about 50 time units the three 
graphs look identical. 

On the other hand, in the bottom part of figured], when we set (3 = 0.5 we 
see that the three oscillators do not synchronize. In fact they seem to tend 
to the type of behavior known as ant i- synchronized or splay state, in which 
each of the three oscillators perform the same motion, but with a phase lag of 
Y' among the oscillators. When many oscillators are involved, such behavior 
will lead to the averaging out of the oscillations, hence no periodicity will be 
observable at the macroscopic level. 

Our aim is to obtain some understanding of phenomena of synchronization 
and desynchronization in oscillators coupled through an environment, such as 
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those demonstrated above. There is an extensive literature concerned with 
the analysis of synchronization of coupled oscillators (see [TTJ [THl EE] and 
references therein), but most studies deal with directly coupled oscillators, 
rather then oscillators indirectly coupled through an environment. In some 
numerical and theoretical studies of systems of the form (11. ip . (11. 21) . 'steady 
state' approximation on (jl.2p is made in order to transform it into a directly 
coupled system: the term y on the left hand side of (11.21) is replaced by 0, the 
resulting algebraic equation is solved for y in terms of x, and the expression 
for y is substituted in (II. ip . Although in some cases this yields a good 
approximation, it is not always so, and it is not hard to find examples in which 
direct simulation shows synchronization of the 'approximate' system, while 
the system (ll.ip . (ll.2p does not synchronize, or vice versa. In this paper we 
study the indirectly coupled system without this steady-state approximation. 

The assumptions that the dynamics of the n oscillators, and their effects on 
the environment, are identical, should be considered as idealizations which 
will only be satisfied in an approximate sense in any real system. However, 
in studying phenomena such as synchronization, it is useful to consider the 
heterogeneity of the characteristics of the oscillators as a perturbation of 
an idealized system of identical oscillators, and the results in this case are 
robust in the sense that results for the idealized systems will also explain 
and predict the behavior of heterogenous systems, at least when the degree of 
heterogeneity is sufficiently small. Indeed, if the system of identical oscillators 
has a stable synchronized oscillation, then by the basic perturbation theory 
of periodic solutions, any sufficiently close heterogeneous system, that is with 
f,g,h in (ll.ll) - (ll.2p replaced by fk,gk,hk, will have a stable periodic solution 
which is 'almost synchronized' in the sense that \xj(t) — Xk{t) \ is small for all 
j,k. 

Another approximation implicit in the model (ll.ip . (ll.2p is that the environ- 
ment is homogeneous, which, in the biochemical context, means that the 
various chemicals diffuse in the environment on a time scale which is faster 
than that of the reactions in the cells and the diffusion across the cell mem- 
branes, or alternatively that the medium surrounding the cells is stirred. 

In section [3] we prove a basic result, Theorem 13.11 which allows to determine 
the (in) stability of synchronized oscillations of (II. ip . (11.21) in terms of the 
stability of two associated linear systems with periodic coefficients. Notably, 
these two linear systems, hence the stability of synchronized oscillations, 
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do not depend on the number n of oscillators. Theorem 13.11 thus reduces 
arbitrarily large problems to a pair of problems which are of fixed size, by 
exploiting symmetry. 

In section H] we make a general study of systems of oscillators which are 
weakly coupled in the sense that they are described by (11. II) . (11. 2p with |/3| 
small. Using Theorem 14 . 1 1 and some perturbation calculations, we prove that 
a system of oscillators coupled through an environment can be synchronized 
using an arbitrarily weak coupling, provided that (3 is chosen to be of the 
'right' sign, and derive a formula which tells us what this right sign is. 

In sections |5] and [6] we present two examples of applications of our analytical 
results to specific systems. 

In section [5] we apply the general Theorem 13.11 to the study of a model of 
periodic release of GnRH, proposed by Khadra and Li |14| . We show that 
whenever this model, for a single cell, produces oscillations, the oscillations 
of any number of cells coupled through the environment will synchronize. 

In section [6] we apply Theorem 14. 11 which deals with the weakly-coupled case, 
to the particular example of indirectly coupled A — uj oscillators, to derive 
explicit conditions for (de)synchronization in the weak coupling regime. 

In section [7] we conclude with some comments on the applicability of the 
results obtained in this paper to the study of systems of the form (ll.ll) - (ll.2p . 

2 Stability, instability and linearization 

We recall the definitions of relevant notions of stability. Let 

z=$(z), (2.1) 

where $ : M. N — > M. N , be a system of differential equations. Let z(t) be a 
T-periodic solution of (12.11) . We denote by O C M. N the corresponding orbit 

O = {z(t) | t G E}. 

z(t) is said to be orbitally asymptotically stable if there exists an open 
set O C U C R , so that for any z$ G U, there exists a p such that the 
solution z(t) of the initial value problem (12. ip . z(0) = z , satisfies 

lim \z(t) -z(t + p)\ = 0. 

t— >oo 
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z(t) is said to be unstable if there exists an open set O C U C R N such 
that for any open set O C V C M N there exists z G V and £ > so that the 
solution z(t) of the initial value problem (I2.ip . z(0) = zq, satisfies z(t) ^ U. 

We now recall the notion of linearized stability. With a periodic solution z(t) 
of (12.11) we associate the T-periodic linearized equation 

w = &(z(t))w. (2.2) 

A Floquet multiplier of (I2.2p is a (generally complex) number /x for which 
(12. 2p has a solution w(t) satisfying 

w(t + T) = iiw(t). 

In other words, defining the fundamental solution of (12.21) as the N x N- 
matrix valued function H(t) satisfying 

H(t) = $'(z(t))H(t), H(0) = I, 

the eigenvalues of H(T) are the Floquet multipliers associated with (12.21) . 
By differentiating (12.11) with respect to t, we get 

I = $'(z(t))l, 

so that z is a T-periodic solution of (12. 2p . which means that fi = 1 is always 
a Floquet multiplier of the linearized equation. The periodic solution z(t) is 
said to be non-degenerate if /i = 1 is a simple Floquet multiplier (that is, it 
is a simple eigenvalue of the matrix H(T)). It is said to be linearly stable 
if it is non-degenerate and all Floquet multipliers other than /i = 1 have 
absolute values strictly less than 1. It will be said to be linearly unstable 
if there is a Floquet multiplier with absolute value strictly greater than 1. 

A fundamental result (see e.g. [20J, Ch. V, Theorem 8.4) states that 

Lemma 2.1 (i) A linearly stable periodic solution is orbitally asymptotically 
stable. 

(ii) A linearly unstable solution is unstable. 

From now on we will refer to linearly (un) stable periodic solutions simply as 
(un)stable. 
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It should be noted that the notions of stability defined above, and the results 
which will be obtained below, are local. Stability of the synchronized oscil- 
lation does not imply that synchronization will be reached from all initial 
conditions, but only that it will be reached with some positive probability - 
that is for a set of initial conditions of positive measure. 



3 Criterion for stability of a synchronized os- 
cillation 

To study stability of a synchronized oscillation Xi(t) = ■■■ = x n (t) = 
x(t),y(t) of 

x k = f(x k ,y), 1 < k < n, (3.1) 

3 " 

y = g(y) + -Y,h{x j , y ), (3.2) 

i=i 

where x(t), y(t) is a T-periodic solution of the single-oscillator system fl 1 . 6 j) . ( TT77T) . 
we linearize the system (13. II) . (13. 21) around this solution, obtaining 

w k = fx(x(t),y(t))w k + f y (x(t),y(t))z, 1 < k < n, (3.3) 

3 " 
z = f -h x (x(t),y(t)) ™i + W(y(t)) + (3h y {x{t),y{t))]z. (3.4) 

n i=i 

Although the system f)3.3p . (13.41) is an nd + p-dimensional one, we will show 
below, using some simple linear algebra, and exploiting the symmetry of 
the system with respect to permutation of the oscillators, that the study of 
its stability reduces to the study of the stability of two linear systems, of 
dimensions d + p and d respectively. When n is large this is a huge reduction 
in the complexity of the problem. 

Theorem 3.1 Let Xi(t) = ■ ■ ■ = x n (t) = x(t),y(t) be a T-periodic synchro- 
nized oscillation of Ii3. ip ./ TOj) . This solution is stable if the following two 
conditions hold: 

(CI) All of the d Floquet multipliers of the T-periodic linear equation 

w = f x (x(t),y(t))w (3.5) 
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have absolute values less than 1. 
(C2) The Floquet multiplier \x = 1 of the T -periodic linear system 
w = f x (x(t),y(t))w + f y (x(t),y(t))z, 

z = {3h x (x(t),y(t))w+[g'(y(t))+f3h y (x(t),y(t))]z (3.6) 

is simple, and all the other d + p — 1 Floquet multipliers have absolute 
values less than 1. 

If one of the Floquet multipliers of either A3. 5\) or A3. 6\) has absolute value 
greater than 1, then the synchronized oscillation is unstable. 

There is an illuminating interpretation of the conditions (C1),(C2). Con- 
dition (C2) says that (x(t),y(t)) is stable as a solution of single-oscillator 
system f ll.6p . fll.7p . As we noted in the introduction, this is a much weaker 
condition then the statement that the synchronized oscillation is stable as 
a solution of the system fl3.ll) . (13. 21) . However Theorem 13.11 tells us that the 
only condition that we have to add in order to get this stronger conclusion 
is (CI), that is the stability of the linear system (I3.5p . This is the system 
one would obtain by looking at x(t) as a periodic solution of the periodically 
forced system 

i = f(x,y(t)), 

with y(t) considered as a given forcing, and asking for the stability of x(t) as 
a solution of this forced system. 

Note that the conditions (C1),(C2) do not depend on n, so that we see that 
stability of the synchronized solution x\ = ... = x n = x, y, where x,y is 
an oscillation of the single-oscillator system fll.6p . fll.7p does not depend on 
n > 2. In other words, if two oscillators synchronize then any number of 
oscillators will synchronize, provided the ratio f3 of the total intracellular 
volume to the volume of the environment is maintained fixed. 

In order to apply Theorem 13.11 to a particular system, one needs to verify 
the conditions (CI) and (C2), so one needs to study the nonautonomous pe- 
riodic systems (I3.5P and (13.61) . a task which may not be easy. Moreover, in 
general one does not have an explicit expression for the periodic oscillation 
x(t), so that even writing down these systems is not possible. Therefore in 
general the verification of these conditions will involve numerical computa- 
tions. There are, however, systems for which the conditions can be verified 
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based on general considerations. Such an example is presented in section 
In addition, Theorem 13.11 is useful for deriving other general results, as we 
demonstrate in the investigation of weakly coupled oscillators in section HI 



Proof of Theorem \3.1\ Written in matrix notation, the system (13.31) . (13.41) is 

/ Wi\ ( Wi\ 

= A(t) 



where 



A(t) 



(a{t) •• 

a(t) 

••• '■ 

■•• 

\c{t) c(t) 



w, 
\ z ) 



(3.7) 



b(t)\ 
b(t) 








a(t) b{t) 
■■ c(t) d(t)J 



a{t) = f x (x{t),y(t)), b(t) = f v (x{t),y{t)), 

c(t) = ^h x (x(t),y(t)), d(t) = g'(y(t)) + f3h y (x(t), y{t)). 

In order to determine the stability of the synchronized oscillation, we wish 
to find the Floquet multipliers of (13. 7J - 

The solution of (13. 1\ given by 



(x\ 

X 



(3.8) 



corresponds to the Floquet multiplier 1. The synchronized oscillation is lin- 
early stable if this is the only solution corresponding to the Floquet multiplier 
1, and the other Floquet multipliers have absolute values less than 1. We now 
display more solutions of (13. 7p and their corresponding Floquet multipliers. 
Assume that w : K —>■ W d satisfies (13.51) . Then by direct inspection one sees 
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that 



/ w \ 
—w 




w 


-w 






/ w \ 



—w 




(3.9) 



V o / V u / v u / 

are linearly independent solutions of (13. 7p . Thus the d Floquet multipliers 
of (13.51) are also Floquet multipliers of (13. 7p . and condition (CI) implies that 
these Floquet multipliers have absolute values less than 1. Since each of these 
Floquet multipliers correspond to n — 1 eigenvectors, we have accounted for 
(n — l)d of the nd + p Floquet multipliers. 

We note also that if w, z is a solution of (13.61) then 



w 

VI 

is a solution of (13.71) . Thus the Floquet multipliers of (I3.6P are also Floquet 
multipliers of (13. 7p . and condition (C2) implies that these Floquet multipliers, 
except the one corresponding to (13.81) . have absolute values less than 1. Since 
(13. 6p is a d + p system we have now another d + p Floquet multipliers. We 
have thus accounted for all nd+p Floquet multipliers of (13.71) . and shown that 
under conditions (CI) and (C2), the system (13. 7p is stable. The argument 
for instability is similar. ■ 



4 Synchronization of oscillators weakly cou- 
pled through an environment 

In this section we consider the system (13. II) . (13. 21) . under the assumption that 
the coupling of the oscillators to the environment is weak, in the sense that 
|/3 1 is small. We will use Theorem I3.1[ and perturbation calculations, to 
obtain information on the (in) stability of synchronized oscillations in this 
regime. 
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Note that, in the biochemical context, in view of (11.31) the weak coupling 
case arises when the volume of the environment is large relative to the total 
intracellular volume, hence the secretion of a synchronizing agent into the 
environment by a cell has only a weak effect on the concentration of this 
synchronizing agent in the environment. 

When (3 — 0, the oscillators do not influence the environment, hence they 
are also uncoupled from each other. We assume that in such the 
environment has a steady state y, and that when the environment is in the 
state y the oscillators have a periodic solution x(t): 

Assumption 4.1 (i) The equation 

y = g{y) (4.1) 

has a stable steady state y, that is g(y) = 0, and that the eigenvalues 
of the matrix 

A = g'(y) 

have negative real parts, 
(ii) The equation 

x = f(x,y) (4.2) 
has a non-degenerate T -periodic solution x(t). 

Part (i) of Assumption 14.11 holds, for example, in the simplest case, in which 
the dynamics in the environment is 'trivial', consisting simply of the decay 
of the various species, so that g(y) = —Dy, where D is a diagonal matrix 
whose diagonal coefficients are the various rates of decay, and y = 0. We 
mention here a paper of Watanabe (see [24J, Theorem 3), which includes 
a perturbation result for synchronized oscillations of systems of the form 
(13. ip . (13. 21) when — > in the case that g = 0, which is excluded by 

Assumption 14.11 

When (3 = the oscillators are uncoupled, so they will not be able to syn- 
chronize, and the most we can expect, in case that x(t) is a stable periodic 
solution of (14. 2p . is that oscillator k will tend to x(t + pu) with different and 
unrelated values of pk- We shall show that for \(3\ > sufficiently small there 
is a synchronized oscillation, but the stability of this oscillation depends on 
the sign of (3: there exists a number a such the synchronized oscillation is 
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stable if <j/3 > 0, and is unstable if a (3 < 0. We shall give an explicit formula 
for computing a. 

Theorem 4.1 Assume that Assumption \4-l\ holds. Then there exists /3q > 
and smooth functions T(P) , x(/3,t), y(j3,t), with 

r :(-#),#)) -(0,oo), 
x: (-#,,/%) xR^M d , y: (-#,,#,) xl^F, 

and 

T(0) = T , x(0,t)=x(t), 2/(0,*) = 2/, 

so £/ia£: 

(%) For a/I |/3| < 

=x 2 (t) = ... = x n (t) =x(0,t), y(t) = y((3,t) (4.3) 

is a synchronized oscillation of the system A3. ( EQ| ) ; wi/i period T(f3). 
(ii) Letting q(t) denote the To-periodic solution of the linear equation 

m = -iW(t),y)Tq(t), (4.4) 

normalized so that 

To 

{q(s),x(s))ds = 1, (4.5) 

we have the following asymptotic formula for the period T(f3) of the above 
synchronized oscillation as (3 — > 0: 

= l-p r [ T °(f y (x( S ),y)e rA h(x(s-r),y),q(t))dsdr + 0(P 2 ). 
J o Jo Jo 

(4.6) 

(^m,) If x(t) is stable as a periodic solution of fiJJty , then defining 

/•To roo 

a= / (fy(x(s),y)e rA h x (x(s - r),y)x(s - r),q(s))drds, (4.7) 



JO J 

we have that, for < |/5| < /3o, ( |^.3| ) stable as a solution of 113. / TOj) «/ 
a/9 > ; and unstable if a (3 < 0. 

(^u) If x(t) unstable as a periodic solution of RJfy , then ft4.3\ ) is unstable as 
a solution of (E2P,(E3) for any < \/3\ < /3 . 
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We recall a fundamental result regarding the perturbation of a non-degenerate 
periodic solution of a differential equation, when the equation is perturbed. 



Lemma 4.1 Let 

x = Q(0,x), (4.8) 

where $ : MxK^ — > M> N is smooth, be a differential equation parameterized by 
the parameter (3. Assume that when (3 = 0, \4-&9 has a To-periodic solution 
x{t), which is non- degenerate. Then there exists 0o > such that if \0\ < 0q 
then h4-8ty has a non- degenerate T \0) -periodic solution x(0,-), where T : 
(— 0o,0o) — > (0, oo) and x : (—0 O ,0 O ) x WL N — > ~R N are smooth functions 
satisfying 

r(o) = r , x(o,t) = x(t), vteM, 

and if x(t) is a linearly (un)stable solution of fl^.flp for (3 = 0, then x((3, ■) is 
a linearly un(stable) solution of fl^.<ff| ) for \(3\ < (3q. 



Proof of Theorem 4-T, The existence of the synchronized oscillation, that is 
a periodic solution of 

x = f(x,y), (4.9) 

y = g(y) + (3h(x,y) (4.10) 

for \(3\ sufficiently small, is an application of Lemma 14.11 Indeed, by As- 
sumption [47TJ when (3 = 0, (I4.9p . fl4.10l) has the non-degenerate T -periodic 
solution x(t) = x(t), y{t) = y, which is perturbed to a T(/9)-periodic solu- 
tion x((3,t),y((3,t) of fl4.9p . fl4.10p for \/3\ small, which gives the synchronized 
oscillation of f!3.ip . 03.21) and proves part (i). 

If x{t) is stable as a periodic solution of (14. 2\ , and since all the eigenvalues 
of A have negative real parts, we have that, for (3 = 0, x(t) = x(t), y(t) = y 
is stable as a solution of fl4.9p . fl4.10p . hence by Lemma l4TTl a;(/3. t),y(3, t) is a 
stable solution of fO|) . flCTD for \0\ < O . Thus condition (C2) of Theorem 
13. II holds for this solution. 

On if x(t) is unstable as a periodic solution of 04.21) . then Lemma I4TT1 implies 
that (C2) does not hold for \0\ < (3 , hence by Theorem 13. II the synchronized 
oscillation is unstable as a solution of f)3.1l) . fl3.2l) . so that we have part (iv) 
of the Theorem. 

To prove part (iii), we need to determine the circumstances under which 
condition (CI) of Theorem 13.11 holds for the synchronized oscillation when 
|/3 1 > is small. 
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For the following computations it will be convenient to normalize the period 
of the periodic solutions to T by setting 

r(P) = (4-11) 

u(/3, t) = x((3, r(/3)t), v(P, t) = y(p, r(p)t), 
so that u(/3,t), v(P,t) are To-periodic with respect to t, and satisfy 

ut(P,t) = T(P)f(u(l3,t),v(0M ( 4 - 12 ) 
v t (P, t) = t(P) \g(v(P, t)) + ph(u{P, t),v(P, *))]• ( 4 -!3) 
To verify condition (CI) of Theorem 13.11 we need, defining 

a(P,t)=T(P)f x (u(P>t)XPM ( 4 - 14 ) 
to check whether all d Floquet multipliers of the T -periodic linear equation 

w = a(P,t)w, (4-15) 

have absolute value less than 1. 

Setting p = in ( 14. 14ft we have a(0, t) = f x (x(t),y), so in this case ( 14.151) 
reduces to 

w = f x (x(t),y)w. (4.16) 

(CI) does not hold for P = 0, since ( I4.16P has the Floquet multiplier \i = 1, 
corresponding to the T -periodic solution w(t) = x(t). However, by our as- 
sumption that x(t) is a stable solution of (14. 2p . all other Floquet multipliers 
of (I4.16P are smaller than 1 in absolute value, and by continuity this remains 
true for (14. 15j) when |/3| > is sufficiently small. Our task, then, is to de- 
termine in what way the Floquet multiplier /i = 1 for P = is perturbed 
when \P\ > is small. Stability of the synchronized oscillation corresponds 
to the perturbed Floquet multiplier being inside the unit disk of the com- 
plex plane. We thus assume that the Floquet multiplier fj, — 1 of (14.151) is 
perturbed to fi(P) for \p\ > 0, where fj,(P) is a real-valued smooth function 
of P with n(0) = 1 - the justification for which is the well-known lemma on 
perturbation of simple eigenvalues (see e.g. [13] . Theorem 5.4). Thus for \P\ 
small there is a solution w(P,t) of 

w t (P,t) =a(P,t)w(P,t), 
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satisfying 

w((3,t + T ) = fi((3)w(p,t), VteR, 
w(0,t) = x(t). 

Denning 

ri(J3)=log(ii(p)), 
w((3,t) = e-^wifrt) 

we have that 

77(0) = 0, 

and w((3, t) satisfies 

w t ((3, t) = [a{/3, t) - v(P)I]w((3, t), (4.17) 

w(P,t + T )=w(P,t), Vt G R, (4.18) 

w{0,t)=x(t). (4.19) 

Stability for small j3 > will hold if r)'(0) < 0, which will imply that r](/3) < 
and hence \n({3)\ < 1. Similarly, stability for small /3 < will hold if rj'(0) > 
0. Part (iii) of the Theorem will be proved by showing that 

?/(0) = -<r, (4.20) 

where a is given by (14. 7p and the rest of the proof is devoted to this compu- 
tation of 7/(0). 

Differentiating (14.171) with respect to f3 we have 

+[a{M - v{P)i\MM 

and putting (3 = and rearranging we get 

w pt (0,t) - f x (x{t),y)w p {0,t) = M0,f)-)/(0)i]i(i). (4.21) 

Taking the inner product of both sides of (14.211) with q(t) (defined as the 
To-periodic solution of (14.41) ) and integrating over [0, To], noting that, using 
integration by parts and (14. 4p we have 

r-T 

(^(0, s) - f x (x(s), y)wp(0, s), q(s))ds 
(Wf,(0, s), q(t) + [f x (x(s),y)}*q(s))ds = 0, (4.22) 
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we get 

°([a p (0,s)-r) , (0)l]x(s),q(s))ds = 0. (4.23) 



o 



Using gSD and fl4T23l) we get 



r/(0) = I °{ap{0,s)x{s),q{s))ds. (4.24) 
Jo 

From (I4.14p we have 

a(P,t)x(t) = r(P)f x (u(P,t),v(P,t))x(t), 
and differentiating this with respect to (3 we get, 

ap(P,t)x(t) = T'{P)f x (u{P,t),v(p,t))x{t) 
+r(P)f xx (u(P, t),v{p, t))[up(j3, t),x{t)\ 
+r((3)f xy W, t),v(j3, t))[Z(t),Vf,(0> *)]• ( 4 - 25 ) 
Putting (3 = in (I4.25P we have 

a^(0,t)^(t) = r'(0)l(t) + f xx (x(t),y)[up(0,t),x(t)] 

+f xy (x(t),y)[x(t),vp(0 1 t)}. (4.26) 

Differentiating (14.121) with respect to (3 we have 

u m ([3,t) = r'(P)x(t)+T(P)f x (u{P,t),v{P,t))u p (P,t) 

+ T(P)f v {u(i3,t)MM)MM, ( 4 -27) 
and setting p = in (14.271) we get 

U/3t {0, t) = r'{0)x{t) + f x (x(t),y)M°> *) + fv(z(t),V)M°> *)■ ( 4 -28) 

Differentiating (14.281) with respect to i we get 

^(0,*) = r'^l^ + Z^^),!/)^),^^,*)] +/ <B (x(*),|/)tt i 9t(0,t) 

+ f xy (x(t), y)[x(t), v p (0, t)\ + fy(x(t),y)vftt(0, t). (4.29) 

Combining ffl~26j) and ffl~29]) we get 

ap(0,t)S{t) = up tt (0,t) - f x (x(t),y)upt(p,t) - f y (x(t),y)vp t (0,t), 
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which we can also write as 

ap(0, t)x(t) + fy(x(t),y)v pt (0, t) = u m (0, t) - f x (x(t),y)up t (0, t). (4.30) 

Taking the inner product of (14.301) with q(t) and integrating over [0,T ], and 
noting that, using (14.41) . the right-hand side vanishes, we get 

To i-Tq 

(ap(0,sfx(s),q(s))ds + / (f y (x(s),y)vp t (0,s),q(s))ds = 0, 

Jo 

which, together with (I4.24p . gives 

V(0) = ~ I °{f y {x{s),y)v m {0, S ),q(s))ds. (4.31) 
Jo 

We now compute vp(0, t). Differentiating (14.131) with respect to j3 and setting 
[3 = 0, we have 

v m (0,t) - Avp{0,t) = h{x(t),y). (4.32) 
Solving (I4.32p for vp(0,t), we get 

vp(0,t) = e tA vp(0,0) + / e rA h{x(t-r),y)dr, (4.33) 



and using the fact that vp(0,t) and x(t) are T -periodic we have 

v p (0,0) = [I-e T ° A )' 1 [ ° e rA h(x(-r),y)dr. (4.34) 



We have 

-To 



/ e rA h(x(-r),y)dr = J2 e {r+kTo)A h(x(kT - r),y)dr 
Jo k=0 Jo 

= [yV ToA ] f ° e rA h(x(-r),y)dr= [J-e^]" 1 f ° e rA h(x(-r),y)dr, 
1 k=o J Jo Jo 

where we have used the periodicity of x(t), and the fact that the geomet- 
ric series converges because the eigenvalues of A have negative real parts. 
Therefore ( I4.34p can be rewritten as 



oo 

rAi 



^(0,0)=/ e rA h(x(-r),y)dr. (4.35) 
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From (I4.33P and (I4.35P we have 

POO pt 

v (O,t) = e tA / e rA h(x(-r),y)dr + / e rA h(x(t - r),y)dr 
Jo Jo 

rt 

e rA h(x(t — r),y)dr + / e rA h(x(t — r),y)dr 

Jo 

e rA h(x(t-r),y)dr, (4.36) 
and, using fl432|) and fl436|) . 



v pt (0, t) = h(x(t),y) + Avp(0, t) (4.37) 

poo 

h(x{t),y) + A e rA h(x(t-r),y)dr 



oo 



A I e rA [h(x(t-r),y) - h(x(t),y)]dr 
d 



oo 

rA 



iirL e rA ][h(x(t-r),y)-h(x(t),y)]dr 



oo 



d 



e rA [h(x(t-r),y)-h(x(t),y)} - e rA -h(x(t - r),y)dr 

r=o J dr 

e rA h x (x(t — r), y)x(t — r)dr. 

Substituting the expression (I4.37P for vp t (0,s) into (14.311) . we get (I4.20p . as 
we wanted. We now prove part (ii). We rewrite (14.281) as 

upt(0,t) - f x {x{t),y)u^,t) = r'{0)x{t) + f y {x{t),y){v p {0, t)). (4.38) 

Taking the inner product of (14.381) with q(t) and integrating over [0, T ] using 
the fact that, due to (I4.4p . the left-hand side vanishes and (14.51) . we get 

r'(0) = ~ [ °(f y (x(s),y)vp(0,s),q(s))ds. (4.39) 



Substituting the expression (14.361) for ^(0, s) into (14.391) . we get 



r'(0) = - / (f y (x(s),y)e rA h(x(s-r),y),q(t))dsdr 
Jo Jo 

which, in view of (14. lip implies (14.61) . ■ 
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5 Application to a model for the pulsatile se- 
cretion of GnRH 



In this section we apply Theorem 13.11 to study a mathematical model, pre- 
sented by Khadra and Li [TJ], whose aim is to explain the synchroniza- 
tion of the periodic (with period approximately 1 hour) secretion of GnRH 
(gonadotropin-releasing hormone) by GnRH neurons in the hypothalamus. 
The explanation proposed for this synchronization phenomenon, based on a 
range of experimental results, is that the GnRH neurons have receptors for 
GnRH, so that the concentration of GnRH in the environment influences their 
dynamics by binding to these receptors, thus inducing an indirect coupling 
leading to synchronization. 

We very briefly describe the mechanisms involved in the model of [13], and 
refer to that paper for details. The binding of GnRH to the receptors on a 
GnRH neuron activates three types of G-protein G s ,G q ,Gi in the cell, and 
the concentrations of their activated subunits a s ,a q ,ai in the cell are de- 
noted by S, Q, I, respectively. a s activates the production of cAMP, whose 
concentration is denoted by A. a q induces the release of Ca 2+ , whose concen- 
tration is denoted by C, from intracellular stores. C and A act in synergy to 
induce the secretion of GnRH from the neuron. G denotes the concentration 
of GnRH in the external environment. These causal relations are modeled 
by the following differential equations: 



where all the parameters, and the nonlinearities Hs, Hq, Hi, F g , Fa, F g , are 
positive. In [T3] these nonlinearities are taken as: 



C = Jin + [I + v c F c {C, Q)){C ER - C) - k c C, 
A = bA + v A F A {S,I) -k A A, 
G = b G + » G F G (C,A)-k G G, 



S = u s H s {G) - k s S, 
Q = vqH q (G) - k Q Q, 
I = uiH s (G) - kil, 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 



H S (G) 



G 4 



Hq(G) 




Kl + G 4 



F C (C,Q) = Q, F A (S,I) 



fuS 



F G {C,A) = {ACf. 



(5.7) 



I + h 



21 



For our results we do not need to assume these specific forms. 

When dealing with n GnRH neurons coupled through the GnRH in their 



environment, the 


model becomes [14) (1 < k < n) 








S k = v s H s {G) - 


ksSk, 




(5.8) 




Q k = vqH q (G) - 


kqQk, 




(5.9) 




i k = ujHsiG) - 


kih. 




(5.10) 


Ck = 


--JiN + [l + vcF c {C k Q k )](C ER - 


Ck) — k c Ck, 


(5.11) 




A k = b A + VAF A {S k i I 


k) ~ kjs^Ak-, 


(5.12) 




G = b G + — Y'F G (C j 
n ' 

3=1 




k G G, 


(5.13) 



where S k ,Qk, h, h,Ck, A k are intracellular concentrations of the various species 
in neuron k, and G is the concentration of GnRH in the intercellular medium. 

We are going to prove results which say that if a single neuron, placed in 
the environment, performs periodic oscillations (in other words if the system 
(I5.ip - fl5.6l) has a stable periodic solution), then a population of such neurons 
will synchronize (in the sense that the synchronized oscillation is stable). 

The following simple lemma will be used. 
Lemma 5.1 Assume a{t) satisfies 

a(t) > a* > 0, Vt > 0, 

and f(t) satisfies 

f(t) + a(t)f(t)=y(t) 

where 

\y(t)\ < me~ kt \/t > 

with k > 0. Then f(t) converges exponentially to as t 
m', k' > so that 

\f(t)\ < m'e- k,t Vt > 0. 



(5.14) 

(5.15) 
oo: there exist 
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Proof: f can be written explicitly as 

f(t) = exp ^ — J a(s)ds S jf(0) + J exp ^ — J a(r)drjy(s)ds. 

Therefore, using fl5.14p . fl5.15l) we have 

|/(*)| < exp(-ta*)|/(0)| +me-"* i / e {a *~ k)s ds. 

Jo 

= exp(— to*) |/(0) | +m- 



^ Jet g CL^t 



which gives the exponential decay. ■ 

We assume that the system fl5.ip - fl5.6l) has a stable periodic solution S(t), 
Q(t), I(t), C(t), A(t), G(t) (that is, we assume that a single neuron performs 
oscillations). Condition (C2) of Theorem 13.11 holds by this assumption. To 
obtain the stability of this solution as a synchronized oscillation of f!5.8p - 
f )5.13p . we need, according to Theorem 13.11 to verify (CI), that is to show 
that the Floquet multipliers of the periodic equation system 

(5.16) 

(5.17) 
(5.18) 



s = 


-k s S 


Q = 


~k Q Q 


I = 


-kj 



c 



I + k c + u c F c (C(t), Q(t)) + -^L(C(t),Q(t))(C ER ~ C(t)) 



C 



+ -^-(C(t),Q(t))(C ER -C(t))Q, (5.19) 

A = VA -^(S(t), I(t))S + VA -^(S(t), I{t))I - k A A. (5.20) 

have absolute values less than 1, or in other words that any solution of this 
system decays to at an exponential rate as t — » oo. Let S(t), Q(t), I(t), 
C(t), A(t), be a solution of fl5TT6|) - fl5T20|) . From fl5T6|) - fl5TT8|) we get that 

§(t) = S(0)e~ kst , Q(t) = Q(0)e- k ^, I(t) = /(0)e~ M , (5.21) 

so that these components certainly decay exponentially. 
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Substituting flEZB into (ET201) we have 

BP ~ BP 

A+k A A = UA -^(S(t), I(t))S(0)e- kst + v A ^f(S(t), /(t))J(0) e -^. (5.22) 

The right-hand side of f 1 5 . 2 2 j) decays exponentially, and k A > 0, hence, by 
Lemma [5. 1[ A(t) decays exponentially. 

Substituting (I5T2TD into fl5TT9l) we get 

-BP 

C + K{t)C = -^-(C(t),Q(t))(C ER - C(t))Q(0)e- k ^, (5.23) 

where 

BP 

K(t) = l + k c + u c F c (C(t),Q(t)) + -g§(C(t), Q(t))(C ER - C(t)). (5.24) 

The right-hand side of ( 15.24ft decays exponentially, so to use Lemma 15.11 in 
order to show that C(t) decays exponentially, we must show that 

mmK(t)>0. (5.25) 



If we assume that Fc does not depend on C, so that 

B P 

^(C,Q) = 0, (5.26) 

which holds in the case of ( 15.71) . then 

K(t) = l + k c + v c F c (C(t), Q(t)) >l + k c >0, 
hence ( 15.251) holds. We thus have 

Proposition 5.1 Assume that 115. 26\) holds and that the system h5. l\) - l[5~b}) 

has a stable periodic solution S(t) , Q(t), I(t), C(t), A(t), G(t). Then, for any 
n > 1, the synchronized oscillation S^it) = S(t), Qk(t) = Q(t), /&(£) = I(t), 
A k (t) = A(t) (I < k < n), G(t) is a stable solution of (TO j) -( fO^) . 

Proposition 15. II covers the case in which the nonlinearities are given by (15. 7p . 
but we now also derive a sufficient condition for synchronization without 
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assuming (I5.26p . As explained in [T4] there is evidence for positive feedback 
of Ca 2+ concentration on Ca 2+ release, so that we assume 



-q§(C,Q)>0. (5.27) 
We note that (15 .4p implies that 

JiN+[l + v C F c (C(t),Q(t))]C ER 



C(t) = => C(t) 



k c + l + u c F c (C(t),Q(t)) 



and since 

n JiN+[l + vcF c (C(t),Q(t))]C ER (Jin + ICer n 

° < k c + l + u c F c (C(t),Q(t)) ~ maX I ko + l ' ° ER 

we have that if C{t) = then 

< C{t) < max( ^±^ ,C fi A (5.28) 
\ k c + 1 / 

so that in particular (15.281) holds at the minimum and maximum points of 
C(t) (recall that C(t) is periodic), hence (I5.28P holds for all t. If we assume 
that 

~, Jin + ICer 
^er > ; —j — , 

then we get C(t) < Cer for all t, hence from (15. 24ft and ( 15.27ft we get that 
(I5.25P holds. We thus obtain 

Proposition 5.2 Assume that j5.21 ) and 

C ER k c > Jin (5.29) 

hold, and that the system ( f5.ip - ( TO]) has a stable periodic solution S(t), Q{i), 
I(t), C(t), A(t), G(t). Then, for any n > 1, the synchronized oscillation 
S k (t) = S(t), Q k (t) = Q(t), I k (t) = I(t), A k (t) = A(t) (1 <k < n), G(t) is 
a stable solution of ( 15.<gj) -( f3T7^) . 



The biologically reasonable values of the parameters given in [14], Cer = 
2.5/xM, kc = 5100 min -1 , Jin = 0.2^, are well within the range satisfying 
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f)5.29p . so that proposition 15.21 applies and assures that synchronization will 
occur, even if the nonlinearity F c depends on both Q and C. 

We remark that Khadra and Li also give a simplified version of their model 
([14j. eqs. 8-9). Proving the stability of synchronized oscillations of this 
simplified model (assuming that a single cell has a stable oscillation), by 
verifying condition (CI) of 13.11 is even easier than in the case of the full 
model, and no restriction on the parameters or nonlinearities is needed in 
this case. 

6 Synchronization of indirectly coupled A — u 
oscillators 

A — uj oscillators are given by the equations 

u = \(Vu 2 + v 2 )u — u(Vu 2 + v 2 )v 

v = oj(\/ u 2 + v 2 )u + A(v u 2 + v 2 )v , (6.1) 

where u), A : [0, oo) — > K are given functions. Introducing polar coordinates 
r, 9 in the u, f-plane, we can write (16. ip as 

f = A(r)r, 9 = u(r). 

It is then seen immediately that if tq > is such that A(ro) = then 

u(t) = r cos(ujQt), v(t) = r sin(ci;ot), 

where u = u;(r ). Moreover, if A'(r ) < then this periodic solution is 
stable, and if A'(r ) > it is unstable. 

In the particular case 

\(r) = l-r 2 , u(r) = 1 + 7(1 - r 2 ) 

we get the Ginzburg-Landau oscillator, which in terms of A = u + iv can be 
written as 

i = (l + z(7+l))A-(l+i7)|A| 2 A (6.2) 

When 7 = this is known as the 'radial isochron clock' [12], or 'Poincare 
oscillator' [6]. 
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We shall assume that 

A(1) = 0, A'(1)<0, 
w(l) = l. 

so that (16.1 1) has the stable periodic solution 

u(t) = cos(t), v(t) = sin(t). 

The existence of an explicitly-known periodic solution facilitates analytical 
study of weakly-coupled A — lo oscillators without resort to numerical calcu- 
lations. A — u) oscillators have been used to address a variety of questions 
related to coupled and forced oscillators - see [25], page 163 for references. 

We will use Theorem 14.11 to study a system of n indirectly coupled A — u> 
oscillators (1 < k < n) 

= \ (yjul + v£)u k -u)(^ju\ + v£jv k + F(u k ,v k )y, (6.3) 

Vk = u(y/ul + vl)u k + \(y/ul + v£)v k , (6.4) 

8 n 

y = -ay + -^2(auj + bvj), (6.5) 

3=1 

where a > 0, in the weakly coupled case \8\ — > 0. We will obtain condi- 
tions on the function F and on the parameters a, a, h which ensure that the 
synchronized oscillation, which exists for \8\ small, is (un)stable. 

The system (jOjl -flOl is of the form (lOi.CT). with d = 2, p = 1, 



x 



f(x,y) 



X(Vu 2 + v 2 )u — u(\/u 2 + v 2 )v + F(u, v)y 
U){y u 2 + v 2 )u + X^u 2 + v 2 )v 

g(y) = -ay 

h(x, y) = au + bv, 
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When (3 = 0, the uncoupled system satisfies Assumption 14.11 the equation 
( 14.1 p has the stable stationary solution y = 0, and the equation (14. 2 j) is the 
A — L) oscillator, with the solution 



cos(t) 
sin(t) 



x(t) = 
We have 

cos(f)[A'(l) cos(<) - c/(l) sin(i)] sin(i)[A'(l) cos(t) - w'(l) sin(i)] - 1 

cos(i)[w'(l)cos(t) + A'(l)sia(t)] + 1 sin(f)[c/(l) cos(i) + A'(l) siu(i)] 



One can check by inspection that 
9(f) 



/ \ J_f A'(l) sin(t) + u/(l) cos(t) 
V «a(*) / 2tt V cos(t) + u/(l) sin(t) y ' 

is the 27r-periodic solution of the equation (14.41) . where the coefficient is 
taken to achieve the normalization (14.51) . 

We compute also 

F(u, v 



fv(x,y) = 
h x (x,y) = ( a b) 

rA u /-/ „\ n\_,-«?f bF(u(s),v(s)) 



f y (x(s),0)e rA h x (x(s-r),0) = e- ay 







f y (x(s),0)e rA h x (x(s - r),0)x{s - r) 

aF(u(s), v(s))u(s — r) + bF(u(s), v(s))v(s — r) 





e -ar 



(f y {x{s), 0)e rA h x {x{s - r), 0)x{s - r), q{s)) 
= e~ ar [aF(u(s),v(s))u(s — r) + bF(u(s),v(s))v(s — r)]qi(s) 
e~ ar 

= [— oF(cos(s), sin(s)) sin(s — r) + bF(cos(s), sin(s)) cos(s — r)] 

x [A'(l)sin(s) W(l)cos(s)], 
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hence 



a 



b r 2n r 

— / F(cos(s),sin(s))[A'(l)sin(s) + cos(s)] / 
2tt Jo Jo 



cos is — re 



- ar drds 



a 

2^ 



2 77 



F(cos(s), sin(s))[A'(l) sin(s) + c/(l) cos(s) 



sinfs — r)e ai 'drds 



1 



a 2 + 1 



(6 - aa)[J(l)C 2 + A'(l)d] + (a + a&)[a/(l)C 3 + A'(1)C 2 ] 



where 



c 2 = — 

2vr 



1 

2^ 

2tt 



2jt 



F(cos(s), sin(s)) sin 2 (s)<is, 



F(cos(s), sin(s)) cos(s) sin(s)<is, 



o 



1 f 2n 

6*3 = — / F(cos(s), sin(s)) cos 2 (s)gLs. 
2tt Jo 



(6.6) 



From Theorem 14.11 we thus obtain 



Proposition 6.1 Assuming A'(l) ^ and a ^ 0, i/iere exists (3q > siic/i 
that for \(3\ < (3q ( fff. 3\) - ji675\) /ias a synchronized oscillation, and we have 

(i) If\'(l) < and 

(6 - aa)[w'(l)C 2 + A'(l)Ca] + (a + a&)[a/(l)C 3 + A'(1)C 2 ] > 0, (6.7) 

where C\,C%,C-z are defined by A6.6\) , then the synchronized oscillation is 
stable for (3 > and unstable for (5 < 0. 

(mJ If X'(l) < and £/ie reverse inequality to ^6. 7\ ) holds, then the synchro- 
nized oscillation is unstable for [3 > and stable for (3 < 0. 

fmj I/A'(l) > t/ien t/ie synchronized oscillation is unstable for all \f3\ < (3q. 

In the case of coupled Ginzburg-Landau oscillators (16. 2p . where we have 
A'(l) = -2, u/(l) = -27, we get 

Corollary 6.1 There exists (3q > snc/i i/iai /or < (3 Q the system of 
coupled Ginzburg-Landau oscillators (1 < k < n): 

ii k = [l- {u\ + v 2 k )]u k - [1 + 7 - -y(u 2 k + v 2 k )]v k + F(u k , v k )y, 
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v k = [1 + 7 - -i{u\ + + [1 - [u\ + t$]i> fcj 

2/ = -ay + - ^Tiaiij + bvj), 

U 3=1 

has a synchronized oscillation, which is stable if 

P[{b - aa)( 7 C 2 + d) + (a + ab)(jC 3 + C 2 )] < 0, (6.8) 

and unstable if the reverse inequality holds, where Ci, 62,63 are defined by 
/ TO) . 

7 Discussion 

We summarize here the basic insights provided by the analytical results given 
by Theorems 13.11 and 14.11 and make some remarks about the possibilities for 
applying these results to the study of specific systems. 

Theorem 13.11 shows that the study of synchronization, that is the determi- 
nation of (in) stability of a synchronized solution of (13.11) . (13.21) . which is a 
system of size nd + p, reduces to the study of the study of two linear systems 
of dimensions d (equation I3.5j) and d + p (equation I3.6p . associated with a 
periodic oscillation of a single oscillator. As we have demonstrated in sec- 
tion [5j for certain systems this result can be used to prove synchronization 
without resort to numerical computations, but in general this will not be the 
case. However, as we have noted, Theorem 13.11 has the following implication 
which is significant in general: if the system (13. II) . (13. 2p with n = 2 has a 
stable synchronized oscillation, then the same is true for any n. If the sys- 
tem with n = 2 is studied by numerical simulation and observed to display 
synchronization, then we are assured that the synchronized oscillation will 
be stable for the system with arbitrarily large n. The caveat must be made 
here that since the notion of stability is a local one, it is possible that for 
n = 2 the synchronized oscillation will be globally stable, but for some larger 
n the synchronized oscillation will only be locally stable. Finding criteria 
for global stability of the synchronized oscillation of system (I3.1I) . (I3.2I) is an 
interesting question for further research. 

Theorem 14.11 provides an understanding of synchronization in the case of 
weak coupling (|/3| small). It is shown that synchronization can occur for 
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arbitrarily weak coupling, provided it is of the 'right' sign, as determined by 
the integral a. It is to be noted that in many cases only a positive value for (5 
makes physical sense, as in the case of coupled cells where (3 is given by (II. 3p . 
in which case the condition for synchronization is a < 0. This criterion can be 
used for a systematic studies of synchronization in the weak-coupling regime, 
in dependence on various parameters. An example of such a study was given 
in section [6] for A — u oscillators. In this case calculations are particularly 
simple, because the periodic solution in the uncoupled case is available in 
closed form. More generally such a study can be performed with the aid 
of numerical computation. Suppose that the coupling function h in (13.21) 
depends on some parameters a = . . . , a m ), and we want to determine 
the subset in the space of parameters a that will lead to synchronization 
for small (3 > 0. We first compute the periodic solution x(t) of (14.21) . for 
example by direct numerical simulation. We then substitute x(t) into the 
formula (14 .7p for a. The dependence of a on a follows from a-dependence of 
h, and the function a (a) can be computed by a numerical integration. The 
surface a(a) = will separate the parameter space into synchronizing and 
non-synchronizing regions in the small parameter case. 
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